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We generalize a discrete-time quantum walk on a line by incorporating a nearest-neighbor inter- 
action into the nonlinear model using the coin operator, which depends on the coin state of the 
previous step. We show that our proposed model has an anomalous slow diffusion characterized by 
the porous medium equation. 
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Cellular automata, discrete models that follow a set 
of rules have been analyzed in various dynamical 
systems in physics, computational models, and theoreti- 
cal biology; well-known examples include crystal growth 
and the Belousov-Zhabotinsky reaction. To simulate 
quantum mechanical phenomena, Feynman Q proposed 
a quantum cellular automaton (the Feynman checker- 
board), and the model, defined generally by Meyer Q, 
is known as the discrete-time quantum walk (DTQW). 
Since the DTQW on the graph is a model of a universal 
quantum computation Q, it is of particular use, espe- 
cially in quantum information [5-7]. Furthermore, the 
DTQW has been demonstrated experimentally in vari- 
ous physical systems (sl-flU . 

As the cellular automaton can be mapped to vari- 
ous differential equations by taking the contentious limit, 
some DTQW models can be mapped to the Dirac equa- 
|, the spatially discretized Schrodinger equa- 
or the Klein-Gordon equation 16]. These 
equations have ballistic transport properties, which are 
reflected mathematically in the one-dimensional (ID) 
DTQW with a time- and spatial-independent coin op- 
erator, i.e., a ID homogeneous DTQW [17]. We con- 
sider here the ID DTQW model. Physically speaking, 
the standard deviation of the homogeneous DTQW is 
a{t) ~ t, whereas the unbiased classical random walk 
has a standard deviation of a{t) ^ \ft. In this Letter, 
we analyze a specific nonlinear DTQW model that incor- 
porates a nearest-neighbor interaction with the coin op- 
erator. We show numerically that the probability distri- 
butions of the nonlinear DTQW model have anomalous 
diffusion characterized by (7^=0.5(0 ~ t^'"^ where is 
the g-generalized standard deviation defined later. These 
dynamics are consistent with the time evolution of the 
self-similar solution [l8| of the porous medium equation 
(PME), which is known to describe well the anomalous 
diffusion of an isotropic gas through a porous medium. 
Furthermore, we show analytically that these dynamics 
do not originate from the interference terms in our model 
whereas the ballistic transport property of the homoge- 
neous DTQW is derived from its interference terms [19|. 



In the homogeneous DTQW, the time evolution of a 
quantum particle (walker) is given by a unitary operator 
U defined on the composite Hilbert space 'Ks^'H.c-, where 
Us '■— span{]j), j € Z} is the walker Hilbert space, and 
"He is the two-dimensional coin Hilbert space. For a uni- 
tary operator U , the quantum state evolves in each time 
step t by 
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where the upper a* (lower 5* ) component corresponds to 
the left (right) coin state at the j-th site at time step t. 
As an example, the time evolution of the DTQW is given 

by 
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The j-th site probability at time step t is given by Pj — 



and 



= 1 is satisfied for each time 



step t. 

As a generalization of Eq. 
DTQW model described by 



([S]). we define a nonlinear 
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with the site-dependent rate function 
9] 
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which incorporates the nearest-neighbor interactions. 
Note that if we set the rate function 5* to g ~ cos 9, 
which is time and site independent, then the model in 
Eq. (jH) reduces to the homogeneous model in Eq. ([3]). 

To study the time evolution of the nonlinear DTQW 
model, the initial state should have nonzero coin states at 
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FIG. 1. The probability distribution of the nonlinear DTQW 
model at t — 10^ step. 



the nearest-neighbor sites. This can be easily understood 
by considering the following example. Let us take (afj, bg) 
as the only non-zero initial state. In this case, the rate is 
(7q = because there is no neighboring state. From the 
map in Eq. (j4]) , we see that the nonzero states a.tt — 1 are 
a^i = — 6q and 6} = a^. This gives g^^ = gl ~ 0, and we 
see that the only nonzero state is (a§, 6q) — {—b\,a^i) = 
(— Cq, — fep) at t = 2. This state at i = 2 only differs in 
sign (or phase) from the initial state. Thus if the initial 
state is concentrated at a single site, then no spreading 
occurs; the state only oscillates around the initial site. 

Figure [T] shows a typical probability distribution of the 
nonlinear DTQW after a long-time evolution. The initial 
state was set as (ag,&g) = (a?,&?) = (1/2, i/2). We note 
that the probability distribution diffuses very slowly and 
does not approach a Gaussian. These features are often 
observed in anomalous diffusion. 

One of the famous anomalous diffusion equations is the 
PME [20], defined by 



d_ 

di 
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where the real parameter m > 1 characterizes the de- 
gree of porosity of the porous medium. It is known that 
the PME can be derived from three physical equations 
for the density p, pressure p, and velocity v of the gas 
flow: the equation of continuity, dp/dt -\- V • (pv) = 0; 
Darcy's law, v oc — Vp; and the equation of state for a 
polytropic gas, p cc p'^ , where 7 is the polytropic expo- 
nent and m = 7 -f- 1. One of the peculiar features of 
the PME is the so-called finite propagation, which im- 
plies the appearance of a free boundary separating the 
positive region (p > 0) from the empty region {p = 0). A 
well-known solution of the PME is the Barenblatt-Pattle 
(BP) solution [3], which is self-similar and its total mass 
is conserved during evolution. The evolutionary behav- 
ior of the BP solution was recently studied in the context 
of generalized entropies and information geometry 21 j. 
The BP solution can also be expressed by Tsallis' one- 
real-parameter (q) generalization of a Gaussian function, 
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FIG. 2. (Color onhne) The probability distribution at t = 10^ 
step running averaged over 10 data sets (light blue line) and a 
fit to the g-Gaussian ((Tjl with q — 0.5 (red dot line) to obtain 
the g-generalized standard deviation (Jqit). 



i.e., the g-Gaussian '22f|. In the case of ID space, the BP 
solution is 



Pq{x,t) (X 
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exp„ 



(7) 



with q = 2 — m. Here, cr^(t) is a positive parameter that 
characterizes the width of the g-Gaussian at time t and is 
similar to the variance cr^^i (t) in a standard Gaussian. In 
other words, the parameter a„(t) characterizes the spread 
of the g-Gaussian distribution [isl [j^: 



agit) cx t3-9. 



(8) 



which reduces to (Tg=i(t) (x ^/t in the limit of q — > 1. 
Note that in the same limit, the g-Gaussian reduces to 
the standard Gaussian, exp (— a;^/iT^^i(t)) , and the PME 
reduces to the standard heat equation dp/dt = d'^p/dx^. 

We performed long-time numerical simulations of the 
nonlinear DTQW model [Eq. (g])] for up to t ~ 10* steps. 
To study the asymptotic behavior, we take running aver- 
ages of the numerical solutions to reduce the influence of 
multiple spikes. The averaged data were fitted with the 
g-Gaussian of Eq. ([7]) to determine the corresponding g- 
generalized standard deviation aq{t), as shown in Fig. [21 
We note that the averaged data at each time step is well 
fitted by the g-Gaussian with g = 0.5. 

The long-time evolution of (Jq{t), plotted in Fig. [31 re- 
veals that the time evolution of the nonlinear DTQW 
model is well characterized by (Tg=o.5(i) ^ ^'''^j which 
is the same time dependency for g = 0.5 of the PME 
[Eq. (HI)]. 

The relationship between our model and the PME 
can be explored u sing the decomposition method of Ro- 
manelli et al. 13, 25l ]. in which the unitary evolution of 
a DTQW model is decomposed into Markovian and in- 
terference terms. We obtain the following map for both 
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FIG. 3. The long-time evolution of the g-generalized standard 
deviation o-q{t) (black dots), which is well fitted by 0-5=0.5(4) ^ 
* (black solid line). 



coin distributions |a*p and 

= (1 - + + 2^i-i5jp/3j, 

(9) 

where the two terms including /3* = 9ie[g*a* (&*)*] are 
interference terms, and 5Hc[z] is the real part of a complex 
number z. 

Neglecting the interference terms and introducing the 
abbreviations — |a*p and i?* = I^P, we get the as- 
sociated Markovian model; 

R*X\ + L'+\ = + (10) 
R]X\ - L]+\ - {2(i*^i + - 1} (i?* - 4)- (11) 

The numerical simulation of the associated Markovian 
model is performed under initial conditions of (i?Q, Lq) — 
{Ri,L\) = 1/4, and the typical probability distribution 
shown in Fig. 2] is well fitted by the g-Gaussian with 
q = 0.0. Furthermore, Fig. [5] shows that the time evolu- 
tion of <Jq{t) of the associated Markovian model is well 
fitted to crq=o.o(i) ^ i" '^^, which again is the same time 
dependency as the PME for q — Q. 

It is known that the classical Markovian model, i.e., 
that without the interference terms of the homogeneous 
DTQW, satisfies the standard heat equation in the con- 
tinuous limit, and consequently the associated asymp- 
totic probability distribution is a standard Gaussian. 
This implies that the ballistic transport property of 
the homogeneous DTQW comes from the interference 
term [l9|. We thus consider the continuous limit |23| 
of the associated Markovian model. 

We introduce the density p{x,t) and current j{x,t) as 

pix,t)=L] + R'^, jix,t) = {R]-^)/Ax, (12) 
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FIG. 4. (Color online) The probability distribution of the 
associated Markovian model Eqs. pU|l and (|lip sd, t = 10^ 
step (light blue line). The fitted parameters of the g-Gaussian 
are q = 0.0 and (7q=o.o(i) = 283 (red dot line). 

0q| 

10^ r - 



10^ T 



10^ 10" 10^ 10" 

t [step] 

FIG. 5. The long-time evolution of the g-generalized stan- 
dard deviation (Jq{t) of the associated Markovian model (black 
dots). It is weU fitted by cr,=o.o(t) ~ t°-^^ (black dot line). 



where Ax is the difference of the nearest-neighbor sites. 
Taking a Taylor expansion of Eq. pU)) , we get 

§iPi-^t) + ^jix,t) + l^pix,t) = 0, (13) 

in the diffusion limit, i.e., the quantity (Ax)"^ / At remains 
constant (set to unity here for simplicity) as At, Ax — 
with the one-step time difference At. In a similar man- 
ner, by expanding Eq. ()11|) and taking the diffusion limit, 
we obtain 

1 d 

2(1 - p{x,tj) ox 

which implies a breakdown in Fick's first law (j cx 
—dp/dx) and is the hallmark of anomalous diffusion. By 
substituting Eq. into Eq. (IT5|) . we obtain the follow- 
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ing nonlinear partial difTerential equation (PDE). 
1 



\^p\^^^)-p\^^^)^pM\- (15) 



Evaluating the asymptotic solution of this nonlinear 
PDE, after a long-time evolution, p{x,t) becomes much 
less than unity. In this long-time limit, we have 1 — p w 1 
and p 



0, and Eq. (jlSp is thus well approximated by 



d_ 

dt 



p{x,t) 



1 



(16) 



which is nothing but the PME in Eq. ([5]) with to = 2 
((7 = 0). We thus conclude that the approximated asymp- 
totic solution of Eq. (IT5|) is a g-Gaussian with q = Q. 
In addition, by applying the asymptotic Lie symmetry 
method [27|, we can show that this result is mathe- 
matically valid [i^l. Therefore, the associated Marko- 
vian model exhibits anomalous diffusion described by the 
PME in Eq. (O with rn^2. 

In summary, we have proposed a nonlinear DTQW 
model Eq. Q in which the coin operator depends on 
the coin states of the nearest-neighbor sites. Our non- 
linear DTQW model asymptotically satisfies the PME 
for TO — 1.5(q = 0.5) and exhibits anomalous slow dif- 
fusion <Jq=a.5{t) ^ i*^'* from the probability distribution 
and the time dependency of the standard deviation de- 
fined in the q-Gaussian distribution. Furthermore, the 
associated Markovian model causes this anomalous dif- 
fusion property. Note that although anomalous diffusion 
was found numerically in another nonlinear model [29] 
and an aperiodic model [30j , the probability distribution 
was not analyzed in those cases. 
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FIG. 6. The power spectrum of the difference between the 
probability distribution of our model and the g-Gaussian with 
q = 0.5. To remove the effects of the expectation value, we 
replace x with x — 36.91 in the g-Gaussian with q = 0.5 [Eq. 
©]■ 



Finally, we would like to comment that there remains 
a unsolved problem concerning on the interference term 



in our nonlinear DTQW model. Consider the difference 
between the probability distribution of our model and 
the (j-Gaussian with q = 0.5, as shown in Fig. [2] the 
power spectrum of this difference exhibits a white noise 
as shown in Fig. [S) This power spectrum divided by the 
physical time scale may remain finite in the asymp- 
totic case, which suggests that our nonlinear model may 
be mapped to the stochastic PME, i.e., the PME plus a 
white noise term, in the continuous limit. This stochas- 
ticity must come from the interference term. The prob- 
lem of extracting the stochasticity from a deterministic 
process has been discussed in another context, that of 
Mori's noise [3l|, and further analysis of this model may 
reveal the origin of the stochasticity. 
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Here, we mathematically show that the probability distribution 
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(1) 



is invariant under an asymptotic Lie symmetry [l[ of the nonlinear partial differential equation Eq. (15) in the main 
text, i.e.. 



dtp 



1 



2(1 -p)2 



dxxP^ - P^dxxP 



(2) 



In Eq. ([TJ, Z{aq=o) = 4crg=o/3 is the normalization factor, and in what follows, the argument of this function is 
omitted where possible and dtp is denoted as pt for simplicity. 

We follow the asymptotic Lie symmetry method and notations in Ref. [l[. Under an infinitesimal transformation 
with the generator 



that is 



X ^ ^dx + rdt + fdp, 



X 1-^ X = X + e ^(x, t, p), 
t ^ t = t + e t{x, t, p), 
p^ p ^ p + eif{x,t,p), 



(3) 



the function p{x,t) is mapped to a new function p{x,t), with 

p{x, t) = p{x, t) + t[ip- pxi ~ ptr] p^p^^^f) ■ 



(4) 



(5) 



By applying this to the probability distribution Eq. ([T]), we see that the transformation X with ^ = ~x leaves Eq. ([T]) 
invariant if and only if 



<P = Px^ 



-PxX 



2x^ 



(6) 



Note that t = rj-t remains unrestricted at this stage because p^''^"^ (x) does not explicitly depend on time t. Conversely, 
the function p(x) is invariant under X = —xdx + Tdt + 2x^ /{Za^^^) dp for any r if and only if p(x) is of the form 
given in Eq. (fT|). 

Following the general procedure for a Lie group analysis of differential equations Q , the second prolongation of X 
is described by 



Y = X + ^tdp, + ^xdp^ + ^xxdp 



(7) 
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The coefScients ^t,"^'^, and "i'xx are defined as follows. Under an infinitesimal transformation of X, the partial 
derivatives are transformed as Px '-^ Px + ^ Pt + £ ^ti and p^x Pxx + £ ^xx- We then readily obtain 

Ax 4 

The coefficients ^P*, and are then obtained by applying the prolongation formula (2.39) in Ref. 0|: 



^i>' = {iPp-Tt)pt = -VPt: (9) 



*'' = 'y5p + (<y5p-^p)pp = ^^^+Px, (10) 
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^"^"^ = (/^xa; + ^ifxpPx + ^ppPx + {^P - '^£,x)Pxx = h Zp^^^;. (11) 

We note that Eq. ^ can be written as 

Pt = Ci {pxf + C2 p,, (12) 

with 

^^=2(1^' ^^^2(W)' (^^^ 
and the asymptotic Lie symmetry condition 

Y{pt - Ci{pxf - C2Pxx) = - 2Cip,*^ - - C'Mpxf ~ C'^VPxx = (14) 

with 

2(i!p)3' ^2 = 9pC2 = ^^^1-^, (15) 

can be written in the following compact form: 

Ao{x,t, p) + Ai{x,t, p)px + A2{x,t, p){px)'^ + A3{x,t, p)pxx = 0. (16) 

When the condition in Eq. is fulfilled, each Ak{k — 0,1,2,3) function must vanish separately in the asymptotic 
limit 

\p{x,t)\^0 for |a:|^cx. (17) 
to imply that the variance crg=o also becomes infinity in the asymptotic limit from Eq. ([T]); 

<Jq=o ^ CO for \x\ CO. (18) 

The function A3 can be expressed as 
which must be nonzero as ct^^o ~^ 00 unless we choose 

ri = -2. (20) 

Making this choice, X becomes 



2x^ 

X ^^dx-2tdt + ^^dp, (21) 



'q=0 

and A3 reduces to 



^3 = n (22) 

(1 - pyZa^g^g 
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Thus, A3 — >• as aq=o 00. 

In a similar manner, Aq,Ai, and A2 are given by 

^°=(l-p)Wo' ^^=(l-p)Wo' ^'=(l-p)Wo' ^''^ 

and all become zero as (Tq=Q — > 00. Therefore, we conclude that the distribution in Eq. ([T]) is an invariant solution for 
the transformation X of Eq. (|2ip. which is an asymptotic symmetry for large |a;| of the nonlinear partial differential 
equation Eq. ([2]). 
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